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qs^ ■ The direct production rates of D-wave charmonia in the decays of Z 

is evaluated. The color-octet production processes Z° ^ 3 D j(cc)qq are 
Ph| shown to have distinctively large branching ratios, the same order of 

magnitude as that of J/ip prodution, as compared with other D-wave 
charmonium production mechanisms. This may suggest a crucial chan- 
nel to test the color-octet mechanism as well as to observe the .D-wave 
charmonium states in Z° decays. In addition, a signal for the 3 Dj char- 
monium as strong as J/ip or ip' with large transverse momentum at the 
c3 ' Tevatron should also be observed. 



PACS number (s) :12.38.5a;, UAO.Gx 



1 



The systematic study of the heavy quark bound systems has played a very important 
role in obtaining information not only on the properties of heavy quarks themselves but also 
on Quantum Chromodynamics(QCD). The property of asymptotic freedom of QCD allows 
one to calculate the production and decays of heavy quark mesons perturbatively at high 
energy scale, while the nonperturbative part can be factored out as the wave functions or 
their derivatives at the origin. 

Recent progress in the this area was stimulated by the experiment results of CDF de- 
tactor at the Fermilab Tevatron. In the 1992-1993 run, the CDF dataQ for the prompt 
production of ip and ip' at large transverse momentum region were observed to be orders of 
magnitude stronger than the lowest order perturbative calculations based on color-singlet 
modelp]], which has ever gained some success in describing the production and decays of 
heavy quarkonia. 

To resolve these discrepancies, Braaten and Yuan suggested|J that parton fragmenta- 
tion represents the dominant source of prompt quarkonium production at high transverse 
monentum (Py > 6 GeVQ), though these processes are formally of high order in the strong 
coupling constant. Over the past few years, the fragmentation functions for S- and P-wave 
quarkonia have been calculated to lowest order0[||]0[|^[Q|9[], and they have subsequently 
been utilized to the issues of quarkonium production phenomenology in hadron collider. 
However, they still underestimate the rates of ip and ip' production more than an order of 
magnitude referring to the observed data at the Tevatron mjTD|. These large discrepancies 
have called in question the simple color-singlet model description for quarkonium and stim- 
ulate ones to seek for new production mechanisms as well as new paradigms for treating 
heavy quark-antiquark bound systems that go beyond the color-singlet model. 

To this end, a factorization formalism has recently been performed by Bodwin, Braaten, 



and Lepage [PH in the context of nonrelativistic quantum chromodynamics(NRQCD), which 
provides a new framework to calculate the inclusive production and decay rates of quarkonia. 
In this approach, the calculations are organized in powers of v, the average velocity of the 
heavy quark(antiquark) in the meson rest frame, and in a s , the strong coupling constant. 

The breakdown of color-singlet model stems from its overlook of the high Fock compo- 
nents contributions to quarkonium production cross sections. The color-octet term in the 
gluon fragmenting to has been considered by Braaten and Fleming[|12|] to explain the 

surplus problem discovered by CDF. Taking < Og ( 3 Si) > and < Og ( 3 Si) > as input 
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parameters, the CDF surplus problems for ip and ip' can be explained as the contributions of 
color-octet terms due to gluon fragmentation. By adjusting the values for matrix elements 
< Of^SJ > and < of ( 3 Si) >, the authors of Refs.|Lj[(Lj got self-consistent values of 
them in the contents of NRQCD. 

Even though the color-octet mechanism has gained some successes in describing the 
production and decays of heavy quark bound systems, especially in explaining the transverse 
momentum spectrum of the quarkonium measured at the Tevatron[|TJ|], it still has a long way 
to go before finally setting its position and role in heavy quarkonium physics. Acturally, by 
now the color-octet matrix elements are extracted only by fitting the data. Moreover, some 



recent studies|TI|[|T6|||T7[ show that there are some inconsistencies about the normalization 
of the color-octet matrix elements in different processes. Therefore, the most urgent task 
among others needs to do is to confirm and identify the color-octet quarkonium signals. 

While the first charmonium state, the J/ip, has been found over twenty years, D-w&ve 
states, given the limited experimented data, have received less attention. However, this situ- 
ation may be changed in both experimental and theoretical investigations. Experimentally, 
there are hopes of observing charmonium .D-wave states in addition to 1 (-^(3770)) in a 
high-statistic exclusive charmonium production experiment ||18|j and bb D-wave states in T 
radiative decays|TP|. Theoretically, the analysis based on NRQCD shows that in D-wave 
quarkonium production the color-octet components play an even more important roles than 
in the S- and P-wave cases. 

Recently, there is some clue for the D-wave 2 charmonium state in E705 300 GeV 
ti^- and proton- Li interaction experiment f20|| . In this experiment there is an abnormal phe- 
nomenon that in the J/ifjir + 7i~ mass spectrum, two peaks at ^(3686) mass and at 3.836 GeV 
(given to be the 2 state) are observed and they have almost the same height. Obviously, 
this situation is difficult to explain based upon the color-singlet model, because the produc- 
tion rate would be proportional to the squared second derivative of the wave function at 
the origin for the D-wave state, which is suppressed by 0(t> 4 ), as compared to the squared 
wave function at the origin for the S"-wave state. However, it might be explained with the 
NRQCD analysis. In NRQCD the Fock state expansion for 3 Dj states is 

\ 3 Dj >= 0(l)\QQ( 3 Dj, 1) > +0{v)\QQ{ 3 P J , 1 8)g > +0(v 2 )\QQ( 3 S 1 , 8 or l)gg >+■■-. (1) 

In fact, considering the suppression due to the derivative of wave function at the origin, 
in production processes, contributions of the three terms written out in Eq.(|T|) have the 



same order in v 2 . So in the quark fragmentation, all of them are of the same order in 
both a s and v 2 . However, in the gluon fragmentation for .D-wave charmonium production 
processes, the S-wave color-octet ( 3 Si,8) production is 0{l/a 2 ) enhanced over the color- 
singlet ( 3 Dj,l) and ( 3 Si,l) production in the short distance perturbative sector because 
the color singlet ( 3 Dj, 1) and ( 3 Si,l) have to couple to at least three gluons. The P-wave 
color-octet process is forbidden by charge conjugation invarance. By this argument, it might 
be easy to understand the E705 experiment data as long as the nonperturbative matrix 
element < 8 D2 ( 3 Si) > is about the same order as < Og ( 3 Si) > in magnitude, which is 
just expected in the framework of NRQCD. 

Of course, this explanation requires the dominance of color-octet gluon fragmentation 
over other production mechanisms. At energies in fixed target experiments like E705, 
the color-octet gluon fragmentation dominance may or may not be the case. Moreover, 
the strong signal of J/ip7c + 7c~ at 3.83QGeV observed by E705 is now questioned by other 



experiments ||21||. Nevertheless, if the E705 result is confirmed (even with a smaller rate, 
say, by a factor of 3, for the signal at 3.83QGeV), the color-octet gluon fragmentation will 
perhaps provide a quite unique explanation for the D-w&ve charmonium production, since 
in all other mechanisms the .D-wave production rate is expected to be much smaller than 
that for the S'-wave states. 

No matter what will be the final situation for the £705 result, it does remind us that 
in the NRQCD approach the production rate of D-w&ve heavy quarkonium states can be 
as large as that of S'-wave states as long as the color-octet gluon production mechanism 
dominantes. This scenario can be tested in many processes, in particular, in the Z° decays. 

Well over 10 6 Z° events have been accumulated at the CERN e + e~ collider at LEP and 
further improvement are expected. This makes it possible to investigate rare decays of Z° 
and to precisely test QCD. Among others the production of charmonium states, in particular 
the 3 Dj states, in the Z° decays will be very interesting in testing the color-octet production 
mechanism. 

Recently a study shows p^] that the leading order color-octet process in ot s , say Z° — > ipg, 
has a relatively small branching ratio because of the large momentum transfer, and this is 
also the case for the D-w&ve charmonium production. The dominant color-octet processes 
for Z° —> 3 Djqq as well as Z° —>■ ipqq begin at order a 2 as shown in Fig.l. Here q represents 
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u, d, s, c or b qurks. From Ref.[p2] we readily have 



T(Z^D jqi ) = T(Z q , f-^ < °> D '™ * {»(! - f ) - 2£ln; 

3o mi 



+ [2Li 2 ( r ^)-2Lz 2 ( r ^) 

- 21n(l + 01ne + 31ne + ln 2 ^](l + 2 } (2) 

X 

in the limit m q = 0, where Li^ix) = — J dt ln(l — t)/t is the Spence function. 

o 

From Eq.dD we can get the branching ratios of Z° —>- 3 Djqq. In the numerical calculation, 



we take p3| [fl3| 



m c = 1.5GeV, Ms Dj w 2m c , a s (2m c ) = 0.253 (3) 

and 

< Ol D2 ( 3 Sx) >^< Of'CSx) >= 4.6 x 10- 3 GeV 3 . (4) 

Here, in NRQCD < Ol D2 {*S x ) > should be of the same order as < 0^( 3 5i) > or 
< Of ( 3 Si) >, and we just take as a tentative value for it, where the value of < Of ( 3 £i) > 



was determined by fitting the CDF data for surplus production of ip' at the TevatronfT3||. 
From approximate heavy-quark spin symmetry relations, we have 

<o^( 3 s 1 ) >^l<ol D *( 3 s 1 )>, 

<Ol D *( 3 S 1 )>* 7 -<Ol D >( 3 S 1 )> . (5) 
Summing over all the quark flavors(g = u, d, s, c, b), one may obtain the deacy widths 

Y^T(Z° -^Dtqq) « 0.7 x lO^GeV, 
q 

J2T(Z° ^ 3 D 2 qq) « 1.2 x 10~*GeV, 
g 

Y, r(^° ^ 3 D 3 qq) « 1.7 x lO^GeV, (6) 

and the fraction ratios 

2.0 x 10 -4 , 



r(z° - 


> 3 Diqq) 


T(Z° 


-> gg) 


r{z° - 


> 3 D 2 qq) 


T(Z° 


-> gg) 


T{Z° - 


+ 3 D 3 qq) 


T(Z° 


-> gg) 



3.4 x 10" 4 , 



(7) 



The dominant color-singlet processes occur as shown in Fig. 2 and Fig.3. Corresponding 
to the quark fragmentation in Fig.2, the branching ratios of 3 Dj production in color-singlet 
processes are 2.3 x 10~ 6 , 3.6 x 10~ 6 , amd 1.7 x 10~ 6 for ,7=1, 2, 3, respectively, which are 



obtained from the universal fragmentation calculations [24]. There should also be color-octet 
processes through quark fragmentation as in Fig.2. However, indirect evidence indicates that 
they are not donimant relative to the color-singlet processes ^M. 

The processes in Fig.3 are more complicated. For in the most important kinematic region 
the virtual gluon is nearly on its massshell, 3 Dj production in the gluon fragmentation color- 
singlet process may be separated to be Z° — ► qqg* with g* —> 3 Djgg. We can estimate the 
partial width following the way in Ref . 1 25| , and the differencial decay rate of Z° — > qqg* may 



then be obtained. With the definition 

F(g* - AX) = n^Pig* - AX), (8) 

the calculation of decay distribution P(g* -* 3 Djgg) for the gluon of virtuality \x is very 
complicated and lengthy (the detailed calculation will be given elsewhere ||26|| ), and in the 
nonrelativistic limit it is proportional to the second derivative of the radial wave function 
at the origin. As in the cases of P-wave charmonium production, g* —^Djgg processes also 
have the infrared divergences involved, which are associated with the soft gluon in the final 
state. Strictly speaking, the divergences can be cancelled in the framework of NRQCD, but 
here we simply deal with it following the way of by imposing a lower cutoff A on the energy 
of the outgoing gluon in the quarkonium rest frame. As discussed in Ref.0, the cutoff A 
can be set to be m c to avoid large logarithms in the divergent terms. 

The decay widths of Z° to color-singlet charmonium state 3 Dj by gluon fragmentation 
can be evaluated via 

Mf 

T(Z°^qqg*;g*^ 3 Djgg) = J dfi 2 T(Z° - qqg*)P(g* ^ 3 Djgg), (9) 
where the cutoff A = m c is transformed into a lower limit on p? min = 12m 2 c . 



In the numerical calculation, taking [23] [27 



a s (2m c ) = 0.253, m c = l.BGeV, |i^(°)| 2 = 0.015GeV 7 , (10) 
and summing over all the flavors q (q = u, d, s, c, b), we obtain 

F(^ - qq) - 4 - 3Xl ° ' 



Among the three triplet states of D-w&ve charmonium, 3 D 2 is the most prominant can- 
didate to discover firstly. Its mass falls in the range of 3.810 ~ 3.840 GeV in the potential 
model calculation [ 28 ] [^J, that is above the DD threshold but below the DD* threshold. 
However the parity conservation forbids it decaying into DD. It, therefore, is a narrow 
resonance. Its main decay modes are expected to be, 



3 D, 



J/^vrvr, 3 /J 2 ^ 3 P j7 (J=1,2), 3 D 2 -+ 3g. (12) 



We can estimate the hadronic transition rate of 3 D2 — > J/ip7T + 7T from the Mark III data for 
^(3770) -> J/^7T + 7r-|30[ and the QCD multipole expansion theory [[H|] |B2| . The Mark III 
data give g r(0(377O) -> J/tpn+n-) = (37±17±8) KeV or (55±23±11) KeV (see also 
Ref. P3). Because the S — D mixing angle for ■0(3770) and -0(3686) is expected to be small 
(say, —10°, see Ref. |3"3"[ for the reasoning), the observed -0(3770) — > J/-07r + 7r~ transition 
should dominantly come from the 3 D\ — > J/ip7r + 7r~ transition, which is also compatible 



with the multipole expansion estimate [32]. Then using the relation |31 

> 3 Si2k) = dT( 3 D 1 -> 3 S 1 2n) 



dT( 3 D 2 



and taking the average value of the T(0(377O) — ► J/^7r + 7r ) from the Mark III data, we 
may have 



T( 3 D 2 -> J/ipir + n- 



T( 3 D 1 -»■ J/ipn + n-) « 46 KeV. 



(13) 



For the El transition 3D 2 — > 3 Pj^f(J = 1,2), using the potential model with relativistic 



effects being considered ]34j], we find 

T( 3 D 2 -> Xdl) = 250 KeV, T( 3 D 2 -> Xc2 7) = 60 #eV, 



(14) 



where the mass of 3 D 2 is set to be 3.84GeV. As for the 3 D 2 



estimate gives p5 



T( 3 D 2 -> 3#) = 12 KeV 
From QT3]), flTJ), and flT5|), we find 

r tot ( 3 /j 2 ) « r( 3 /j 2 ^ j/Vw7r)+r( 3 D 2 -^Xci7)+r( 3 D 2 

« 390 KeV 



3g annihilation decay, an 

(15) 



Xc27 ) + r( 3 i} 2 ^3<7) 



(16) 



and 

B( 3 D 2 -> J/tpn+n-) w 0.12. (17) 

Compared with -B^' — * J/ipir + ir—) = 0.324±0.026, the branching ratio of 3 D 2 — > J/if>7r + 7i~ 
is only smaller by a factor of 3, and therefore the decay mode of 3 D 2 — > J/ip7T + 7T~ is 
observable, if the production rate of 3 D 2 is of the same order as if)'. 

The other two states of the triplet, 3 Di and 3 D 3 , are above the open channels threshold 
and are not narrow, and therefore are difficult to detect. It might be interesting to note that 
the OPAL Collaboration at LEP has analysed the J /ipir + ir~ spectrum in Z° decays |BS], and 
there seems to be some events above the background around ?>.77GeV , whether these events 
are associated with the D-w&ve 1 charmonium state ^(3770) might remain interesting. 

In conclusion, from the calculations and discussions above one can clearly see that the 
branching ratios of gluon fragmenting to color-octet 3 Dj states are two to three orders 
larger than the dominant color-singlet processes. This large divergences are much helpful in 
distinguishing the two production mechanisms in experiment. On the other hand, because 
the production rates of if)' and 3 D 2 (cc) in color-octet mechanism are of the same amount of 
magnitude, the 2 charmonium production in Z° decay may provide a crucial channel to test 
color-octet production mechanism at LEP with present luminosity. As for the bb system,, 
the production of D-w&ve bottomonium states by color-singlet and color-octet mechanisms 
are similar to that of charmonium, but the disparity of the branching ratios between these 
two mechanisms is not as evident as that of charmonium production. 

Finally, we would like to point out that the above discussion also applies to the 3 D j 
production at the Tevatron. As in the case of J/if> and if>' production at large momentum 
transverse, the 3 Dj production will also be dominanted by the color-octet gluon fragmenta- 
tion. This implies that a signal for 3 Dj states, which can be as strong as J/if> or if)', should 
be observed at the Tevatron. This will also be a crucial test of the color-octet mechanism. 
Detailed analysis will be given elsewhere. 
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Figure Captions 



Fig.l. One of the contributing Feynman diagrams of color-octet mechanism in Z G ^ 3 Djqq 
processes. 

Fig. 2. One of the Feynman diagrams corresponding to quark fragmentation processes in Z° 
decays. 

Fig. 3. Diagrams for 3 Dj production from gluon jet in color-singlet mechanisms, (a)virtual 
gluon production in Z° decays (b) 3 Dj production in gluon fragmentaion. 
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